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In addition to the final-state interactions, J/ip nuclear production might also be affected 
by parton energy loss. Using the upper limits from Drell-Yan data at SPS and Fermilab 
energies, we estimated energy loss contribution to J/ip production in p-A collisions. The 
results indicated that the effects might be sizeable at 200 GeV while remaining small at 
higher energies. 



1 Introduction 

Charmonium production is one of the cleanest hard probes studied in high-energy hadron- 
hadron, hadron- nucleus, and nucleus-nucleus scatterings. If we look at the reaction in the center 
of mass frame, the parton model tells us that the production results from the fusion of a beam 
parton with energy-momentum fraction x\ and a target parton with energy-momentum fraction 
X2 which are converted into a cc pair that eventually gives the observed charmonium. In the 
absence of nuclear effects the integrated production rate in A\-A<i collision would be given 
by A\ x A2 times that in p-p reaction. The deviation from this is often interpreted within a 
conventionnal absorption model. 

Recently the E866 collaboration provided the first evidence for a distinction between J/ip 
and ip' in the differential production ratioli 

da(p + W ^iP + X)/dx F m 
Ada(p + Be^ ip + X)/dx F 1 ' 

The distinction occurs in the kinematical region where the charmonia have relatively small 
energies E c5 in the target rest frame. This is a definite departure from the observation of 
identical ratios for larger cc energies. It is tempting to describe this result by saying that for 
small cc energies the nucleus "acts as a detector" for the formation of the charmonium states. 
A possible interpretation of this is that for large cc energies a common precursor to J/tp and ip' 
passes through the target and final-state formation takes place outside the nuclear volume. At 
low cc energies formation occurs within the target, leading to a weaker absorption for J/ip than 
ip 1 due to the smaller size of the former. At the time E866 data became available we decided to 
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take these ideas at face value and to test them in a simple but quantitative scenario of nuclear 
absorption I. 

A direct consequence of such an approach is an X2 scaling. At high energy the charmonium 
energy amounts to E c5 = X\E in the target rest frame where E is the beam energy in this frame. 
In terms of the charmonium invariant mass M c5 one has E cS = M^j (2m v X2) . Consequently, if 
the production ratio is in actual facts driven by E c ^ one expects R to be independent of the 
beam energy at a given X2 (3:2 scaling). 

Unlike what we expected, a comparison of J/ip differential production between 800 GeV 
(E866)0 and earlier less precise 200 GeV data (NA3)i as a function of X2 seems to indicate that 
cc production does not follow an xi scaling. However, it could happen that NA3 data reveal 
that at intermediate energies there is a mechanism involved in cc nuclear production beside the 
above mentionned nuclear suppression. 

In fact, it has been argued for a long time that an incoming parton might lose some energy 
while scattering through nuclear matter. Its momentum fraction is then shifted from xi+Axi 
to x\ at the point of fusion resulting in a modification of the effective beam parton density 
by a factor F{x\ + Ax\)jF{x\). Because the parton distributions F{x) drop dramatically at 
intermediate and large x, the effective density of beam partons and consequently the production 
ratio might decrease significantly. 

The aim of these proceedings is to investigate whether energy loss could be a relevant mech- 
anism to explain a breakdown of X2 scaling in J/ip production. 



2 Methods 



Since a corresponding quark energy loss may be expected in Drell-Yan (DY) pair production 
which is moreover not overwhelmed by final-state interactions, we first used the DY process to 
extract limits on parton energy loss. Energy loss is then combined with cc absorption to examine 
the consequences in the J/ip channel. 

To be more specific we chose to test two alternative energy loss scenarios. The first is that 
put forward by Gavin and Milana (GM)S. It predicts a beam parton shift in momentum fraction 

Axi = KiXiA 1/3 , (2) 

where K\ is a free parameter |] We stress that the shift depends on kinematics through x\. 

The second model was proposed by Brodsky and HoyerH in order to fix a potential defect 
in GM. On general grounds they showed that 

Axi < ?l A i/3 ) ( 3 ) 

with K2 ~ 0.5 GeV 2 , an upper bound with which Ax± as given by Eq. (|J) is incompatible 
at large energies. For accessible energies it also sets an upper bound on K\. In relation with 
these considerations an alternative scenario, denoted BH in the following, is then obtained by 
assuming a shift in x% that saturates the bound given by (||), i.e., 

Axi = m A ifl 

s 

where it should be noted here that the dependence on kinematics occurs via the center of mass 
energy. 

*Ki might depend on Q 2 and we needed its value in the region around m^. It turns out that the DY data 
are not precise enough and do not probe a sufficient range in Q 2 to make any definite statement about the latter 
dependence. We thus simply disregarded it. 



The free parameter in each model was fixed with DY data. The data set encompassed, on 
the one hand, the 150 and 280 GeV pion beam data on A = 1 and 195 targets! and, on the 
other hand, the 800 GeV proton beam data on A = 9 and 1840. At 800 GeV the data set spans 
a range in x 2 where shadowing is known to be sizeable. We thus used data points corrected for 
shadowing as given in Ref. 7. At SPS energies the data points lie in the (weak) antishadowing 
region and the corresponding correction is small. Calculations were done using MRST parton 
distributions Is and shadowing corrections were given by the EKS parameterization a. 

The models were then used to assess the possible importance of energy loss in J/ip production 
at 200 and 800 GeV. We assumed that J /if) production results from gluon fusion and took into 
account the color factor correction to energy loss : Ax g = \Ax q [[ Subsequent energy loss was 
combined with nuclear absorption extracted from 800 GeV measurements and the result was 
compared to 200 GeV data!. In the J/ip sector shadowing corrections were not considered (see 
discussion at the end of Section 3). 

3 Results and Discussion 

Drell- Yan production 

In view to quantify parton energy loss, we first considered Drell- Yan dimuon production. At 
800 GeV the wide kinematic acceptance (0.2 < x\ < 0.95) and high statistics of E866/NuSea 
data make it possible to put stringent constraints on parton energy lossQ. The fitted parameter 
values were 

Kl = 3 x 10" 4 (x 2 / n df = 0.69), 

for GM and 

k 2 = 0.08 GeV 2 (x 2 / ndf = 0.72), 
for BH. Their la upper limits, i.e., x 2 / n df= 1, were 

Ki + Aki = 11 x 10~ 4 , 

and 

k 2 + Ak 2 = 0.77 GeV 2 . 

This suggested that the Drell- Yan data — corrected for shadowing — are consistent with zero 
energy loss within both models, but the la upper-limits imply that energy loss effects at a level 
of a few hundreds of MeV/fm cannot be excluded. 

Similarly, limits were estimated from a fit to the ir~ beam NA3 data at 150 and 280 Gevi 
(Table |). The lowest x 2 /ndf were reached with essentially no energy loss for both models 
- with or without inclusion of shadowing. Table p] also indicates that the consideration of 
shadowing effects tended to improve the fits and correlatively to raise the upper limits of parton 
energy loss by a factor of four. 





GM loss 


BH loss 


hi (xl0~ 5 ) 


x7 ndf 


k 2 (GeV^) 


xVndf 


no shadowing 


0.0 (1.0) 


0.86 


0.0 (0.2) 


0.86 


with shadowing 


0.2 (3.7) 


0.43 


0.0 (0.9) 


0.43 



Table 1: Energy loss parameters Ki and Ki from a fit to NA3 data, with and without shadowing corrections. The 
la upper limits ki + Am and K2 + Aac2 are given between brackets. 



'At large zi, qq annihilation into cc becomes predominant t_l. We checked that taking this channel into account 
did not change the overall results. 



Fitting the nuclear dependence of Drell-Yan production in p - A and tt - A reactions leads 
to the conclusion that quark energy loss has to be small if any. Contrary to what we expected, 
the comparison between 200 and 800 GeV data did not allow us to discriminate between the 
different energy behaviours of GM and BH scenarios. A closer look at the data showed that this 
is mostly due to the large error bars in NA3 data. 

J/ip production 

In the parton model, production of cc pairs proceeds via parton fusion and might therefore 
be affected by initial-state effects (e.g., energy loss and shadowing), as well as final-state inter- 
actions. From the preliminary Drell-Yan study, we assumed K\ to be equal to 0.001 and k<i to 
be equal to 0.5 GeV 2 to evaluate the possible energy loss contribution to J/ip production. 
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Figure 1: Energy loss effects on J/ip production according to GM loss (dotted) and BH loss (dashed) at 200 GeV 
(left) and 800 GeV (right). The production ratio is defined as R(A/B) = Ba(A)/Aa(B). 

Figure | shows the J/ip production ratios R(Pt/p) at 200 GeV and R(W/Be) at 800 GeV. 
The most striking result is the large suppression seen at SPS according to BH mechanism whereas 
this suppression remains negligible at Fermilab energies. The J/ip production at 200 GeV is 
already affected (R ~ 0.9) at low x\, and decreases when x\ increases. The situation is different 
in the GM approach where the loss scales with x±, independently of the beam energy. This can 
be seen in Figure [l]j] which also demonstrates that the effect is moderate for k\ = 0.001 with 
a decrease from 1 at small x\ to 0.8 at x\ = 0.7. The decrease seen in both approaches when 
going from small to large x\ is the behaviour that might be reflected by NA3 J/ip data. 

Once energy loss in the J/ip sector was examined, J/ip suppression was calculated taking in 
addition nuclear absorption into account. The results are compared with experimental data! 
on Figure ||| as a function of xp = x\ — X2- First, Figure § reminds us that the model de- 
scribed in Ref. 2 for charmonium suppression (solid curve) is unable to describe large xp NA3 
data. Second, the energy loss mechanism proved to reduce the disagreement with large xp data 
points. In particular, GM loss did not affect low xp production ratios and somehow decreased 
J/ip production at xp ~ 0.5. By contrast, the BH scenario could almost reproduce the large 
suppression seen at large xp but tended to give too strong a suppression at small xp. Though 
it is true that a first level comparison of both models would favor the GM approach, we should 

*The difference between the two dotted curves in Figure |l] arises from the normalization with p and Be targets 
at 200 GeV and 800 GeV respectively. 
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Figure 2: Ratio of J/ip's produced in p(200 GeV)+Pt and p(200 GeV)+p versus xf from NA3 as compared to 
theoretical calculations of the nuclear suppression (solid), with additional GM energy loss (dotted) or additional 

BH energy loss (dashed). 



however stress that the results depend strongly on the values chosen for the parameters k\ and 
K2- This prompted us not to draw too quantitative conclusions. 

In summary, the results obtained may indicate that NA3 J/tp data result from a possible 
interplay between energy loss and nuclear suppression. In particular, the general trend of data 
is better reproduced when parton energy loss is taken into account. Consequently, the lack of 
X2 scaling exhibited by NA3 and E866/Nusea data could be — at least partly — explained with 
such a mechanism. One limitation nevertheless arises from the remaining discrepancy between 
data and theoretical calculations. Further to that, we should mention that other mechanisms 
have been neglected in this study. In particular, gluon antishadowing calculated with EKS98 
leads to an enhancement of J/ip production in platinum near x<i ~ 0.1 (i.e. xp ~ 0.2 for NA3). 

It is hardly necessary to repeat that precise measurements of both Drell-Yan and charm 
differential production over a large kinematic window is mandatory to disentangle all the mech- 
anisms involved in J ftp production. 
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